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Introduction 

Variable  Speed  Control  Moment  Gyroscopes  (VSCMGs),  defined  by  Schaub  in  [1],  effectively 
applied  to  combined  energy  storage  and  attitude  control  in  [2, 3],  and  further  studied  by  Yoon  and 
Tsiotras  in  [4,5],  can  avoid  attitude  singularities  by  varying  wheel  speed.  However,  this  device’s 
key  advantage  lies  in  torque  amplification  via  gimballing  (i.e.  operating  in  its  CMG  mode).  The 
primary  limitation  of  this  CMG  mode,  though,  is  the  existence  of  gimbal-lock  singularity  states. 
Avoiding  escapable  singularity  states,  studied  extensively  by  Wie  for  CMGs  [6]  and  Schaub  for 
VSCMGs  [7],  is  exacerbated  in  the  case  of  simultaneous  attitude  and  power  tracking  require¬ 
ments  [5].  Although  [5]  gives  an  algorithm  for  combined  energy  storage  and  attitude  control  with 
VSCMG  singularity  avoidance,  it  does  not  directly  evaluate  this  algorithm  with  respect  to  other 
robust  CMG  singularity  avoidance  methods  from  the  literature.  For  this  reason,  two  practical  steer¬ 
ing  laws  were  adapted  from  previous  results  for  this  challenging  case  and  compared  to  the  results 
given  in  [5],  First,  Wie’s  Generalized  Singularity  Robust  (GSR)  steering  law,  originally  developed 
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for  attitude  tracking  with  CMGs  [8]  was  modified  to  use  VSCMGs.  Although  the  GSR  effec¬ 
tively  perturbs  gimbals  through  locked  configurations  when  compared  to  classical  CMG  methods 
(see  [6, 9, 10]),  it  also  imparts  small  attitude  torque  errors  which  can  linger  for  a  few  minutes  af¬ 
ter  a  CMG  singularity.  This  problem  can  be  disruptive  to  some  space  missions  that  use  satellite 
off-pointing  (e.g.  precision  earth  imaging)  where  accurate  trajectory  tracking  during  operations  is 
paramount.  However,  VSCMG  wheel  speed  freedom  permits  the  modified  GSR  to  remove  these 
errors,  but  is  limited  in  the  power  tracking  case.  Further  extending  this  modified  GSR  to  track 
power  in  the  null  subspace  of  the  VSCMG  steering  dynamics  (as  in  [2, 3])  enables  an  improved 
steering  law  that  overcomes  CMG  singularity,  eliminates  small  attitude  errors,  and  tracks  a  desired 
flywheel  energy  storage/drain  profile.  Second,  Schaub’s  VSCMG  with  null  motion  algorithm  [7], 
designed  to  exploit  its  CMG  mode  while  more  efficiently  accommodating  gimbal-lock  singularity 
than  previous  VSCMG  approaches,  was  modified  for  the  power  tracking  case.  As  demonstrated 
here,  the  adapted  GSR  and  the  modified  VSCMG  plus  null  motion  algorithms  perform  comparably 
(showing  a  small  improvement)  over  Yoon/Tsiotras’  results  from  [5]. 


Simultaneous  Attitude  and  Power  Tracking  Model 

As  shown  in  [3],  a  model  encompassing  the  spacecraft  dynamics,  kinematics,  required  attitude 
reference  torque,  and  VSCMG  attitude  plus  power  steering  is 


Nd  —  Isc6)-\-  d)/jC(0-t-  B8  -|-  Ef>l  +  -f-  F£l 
Nr  =  klsctf  (P)Pr4~  7f/sC(Ct)  COr)  G)/V(00 


m  = 

Bl  +  ESl  +  Dh  +  FQ.  = 
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(3) 

(4) 
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where  Nd  is  the  inertial  torque  acting  on  the  combined  system  of  the  spacecraft  platform  plus  the 
VSCMGs,  l sc  is  the  system  inertia  assuming  the  gimballed  VSCMG  inertia  (from  rotor  dynamics) 
is  negligible,  P  is  the  Euler  parameter  set  (po,Pi,P2>P3)  representing  the  satellite  body  frame 
orientation  with  respect  to  the  inertial  frame,  CO  is  the  angular  velocity  of  the  body  frame  with 
respect  to  the  inertial  frame  expressed  in  body  frame  coordinates,  (3  and  d)  are  the  body  frame 
time  derivatives  of  P  and  to,  pr,  and  cor  are  command  reference  versions  of  P  and  co,  8  is  an  n  x  1 
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column  matrix  of  n-VSCMG  gimbal  angles,  8  and  8  are  n  x  1  gimbal  angular  velocity  and  angular 
acceleration  column  matrices,  d)  is  a  3  x  3  skew-symmetric  matrix  using  the  elements  of  (0;  B,  D, 
E,  F  are  3  x  n  matrices  transforming  actuator  parameters  from  the  gimbal  coordinate  frame  to  the 
body  frame  where  B  is  a  function  of  8  and  is  constant  in  the  body  frame,  D  depends  on  CD,  8,  and 
O,  varies  in  the  body  frame,  and  is  approximated  as  its  wheel  angular  momenta  component,  E 
transfers  the  wheel  acceleration  to  the  body  frame  and  is  dependent  upon  the  gimbal  angles,  and 
F  depends  on  the  gimbal  angles  and  body  angular  velocity  co  and  is  equivalent  to  the  G)h  defined 
in  [5,6],  lWSd  is  a  MX  «  diagonal  matrix  of  VSCMG  wheel  spin-axis  inertias.  Pr  is  the  required 
power  to  store/drain  in  the  wheels,  while  Pw  is  the  actual  power  contained  in  the  wheels.  Since 
V  =  ClTIWSd,  then  Pw  =  VCl.  Finally,  k  is  a  positive  gain  scalar,  K  is  a  3  x  3positive  definite  gain 
matrix,  and  Nr,  which  uses  k  and  K,  is  the  required  torque  for  stable  spacecraft  attitude  tracking. 

Noted  often  in  the  literature  (e.g.  [1,3]),  CMGs  typically  exploit  torque  amplification  properties 
through  gimbal  rate  control  (i.e  keeping  gimbal  motor  torques  small  such  that  8  «  0).  The  result  is 
a  velocity-based  steering  law  to  replace  eq.  5  (assuming  no  external  torques  are  applied) 


£0  +  Z)8  =  Nr-FCl 


This  velocity-based  steering  law  permits  different  combinations  of  gimbal  rate  and  wheel  accel¬ 
eration  (similarly  torque)  to  achieve  a  desired  attitude  maneuver.  One  can  store  (drain)  energy  by 
increasing  (decreasing)  the  wheel  speeds  in  combinations  that  still  permit  the  desired  net  torque 
for  attitude  tracking.  Eq.  6  captures  the  goal  for  this  additional  requirement:  get  the  actual  power 
in  the  wheels,  Pw  =  VO  to  match  a  time-varying  power  profile,  i.e.  Pw  =  Pr{t).  From  eqs.  6  and  7, 
one  can  define  the  set 
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where  u  = 
and  S=[V 
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,  Nc  =  Nr-  FCl,  Pc  =  Pr{t),  Q—[E  D),  0\xn  is  a  1  x  n  matrix  of  zeros, 


Oix/v]-  Eq.  8  enables  simultaneously  tracking  a  stable  attitude  reference  and  a  power 


Attitude  Singularity  Avoidance 

The  fundamental  VSCMG  attitude  steering  problem  is  ensuring 
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which  yields  a  general  VSCMG  steering  law 

*  =Q'NC  (10) 

o 

L  J  c 

where  Q 1  is  some  generalized  inverse  of  Q  such  as 

=  WlQT  (QWiQ7)-'  (11) 

In  eq.  1 1 ,  a  weighting  scheme  is  used  to  weight  the  gimbal  and  wheel  command  modes,  thereby 
permitting  a  method  for  the  wheels  to  take  control  near  gimbal  singularity  and  gimbal  torque 
advantages  to  prevail  everywhere  else.  One  way  to  combat  this  is  to  use  the  Generalized  Singularity 
Robust  (GSR)  steering  law  introduced  by  Wie  et.  al.  in  [8, 1 1] 

4  =Q*NC  02) 

o 

L  J  c 

where 

Q*  =  WxQt  (QWxQt  +  G)~l  (13) 

This  method  is  adept  at  rescuing  the  system  from  gimbal  lock,  but  does  impart  residual  attitude 
torque  error  during  the  maneuver  and  for  some  time  afterward.  These  errors  can  disrupt  the  ability 
to  track  a  desired  attitude  during  precision  imaging  operations.  One  solution  is  to  use  wheel  control 
to  eliminate  these  errors.  This  is  done  as  follows.  Let  utrue  be  the  commanded  u  required  to  achieve 
Nc,  ugsr  be  the  u  associated  with  the  GSR,  and  uerror  be  the  commands  necessary  to  eliminate  the 
tracking  error  caused  by  the  GSR.  Then  the  following  equations  apply 


Utrue 

—  ugsr  uerror 

(14) 

UGSR 

II 

(15) 

QUerror 

error 

(16) 

Uerror 

=  tfNcenor 

(17) 

Define  Q*  =  WzQT  {QWjQ7)  \  where  the  weighting  matrix  W2  permits  weighting  between  the 
wheels  and  gimbals  independent  of  the  weights  in  W\  and  only  acts  during  singularity.  Next,  derive 
the  error  torque,  NCerror,  where  I3  is  the  3  x  3  identity  matrix,  as 

Ncerror  =  Nc,rue  ~  lVCcsfl  =  Nc  -  QQ*NC  =  (/3  -  QQ*)  Nc  (18) 
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Finally,  assemble  the  components  of  utrue,  which  now  becomes  the  GSR  correction  steering  law, 
UGSRcorr  for  VSCMGs 


UGSRcorr  =  "'rue  =  Q*NC  +  Q*  (/3  -  QQ*)NC  =  Q°NC 


(19) 


Thus,  Q°  =  Q*  +  Qt  (h  ~  QQ*)-  Clearly,  a  logical  choice  for  W2  (where  In  is  the  n  x  n  identity 
matrix)  is 


W2  = 


In  On 
0N  W?IN 


(20) 


—/o  del  I  - 

where  IV2  =  \imdiv>sdj  jt  permits  de-weighting  the  gimbal  mode  when  constmcting 

the  attitude  error  correction  torque  (the  gimbals  are  in  singularity  when  W2  acts).  This  thereby 
promotes  the  wheel  mode  as  desired.  In  contrast,  the  choice  for  Wj  comes  from  typical  applications 
such  as  found  in  [5] 


Wi 


In  On 
On  W* 


(21) 


where  Wg  =  Wgoe  Although  this  VFi  matrix  will  work  in  theory,  it  is  not  practically  efficient 

at  keeping  wheel  and  gimbal  accelerations  small  during  singularity  avoidance. 


Singularity  Avoidance  with  Power  Tracking 


The  technique  to  project  power  in  the  nullspace  of  the  dynamics  matrix  (Q  in  this  case)  was 
shown  in  [3].  The  steps  are  given  in  the  Appendix.  This  yields  the  extended  GSR  simultaneous 
attitude  and  power  tracking  steering  law,  unew, 


Unew  —  ^GSRcorr  "b  —  Q  He  ~b  I*nS  nS  )  (Pc  SQ  Ac) 


(22) 


3.0  Numerical  Example 

The  utility  of  the  extended  GSR  (eq.  22)  and  the  modified  VSCMG  with  null  motion  algorithms 
for  attitude  and  power  tracking  are  best  viewed  in  the  context  of  a  practical  example  such  as 
that  from  [5].  In  this  scenario,  attitude  and  power  tracking  for  five  cases  are  investigated,  I:  no 
singularity  robust  steering,  II:  GSR  steering  with  no  attitude  error  correction.  III:  extended  GSR 
with  attitude  error  correction,  IV:  modified  VSCMG+null  motion  from  [7],  and  V:  VSCMG  attitude 
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and  power  tracking  from  [5].  The  parameters  used  for  these  cases  are  listed  in  Table  1  while  Table 
2  highlights  some  performance  results.  More  specifically,  it  compares  attitude  error,  momentum 
wheel  mode  torque  generated  by  each  VSCMG  (NSmax),  top  wheel  speed  (in  RPM),  nominal  wheel 
speed,  and  the  percentage  difference  between  the  maximum  and  nominal  wheel  speeds.  These  are 
maximum  values  of  respective  physical  quantities  during  the  maneuver.  However,  attitude  errors 
are  compared  only  during  the  tracking  phase  and  the  initial  transient  phase  is  ignored.  One  can 
see  the  attitude  error  is  increased  in  case  II  which  is  then  corrected  in  cases  III,  IV,  and  V.  In  the 
overall  scenario,  smaller  values  of  torque  are  required  for  cases  II  through  V  as  compared  to  case  I 
as  well  as  smaller  variation  in  wheel  speed  for  the  same  cases.  These  results  illustrate  the  improved 
efficiency  of  the  proposed  singularity  avoidance  plus  attitude/power  tracking  methods  (cases  III, 
IV,  and  V). 

The  ensuing  figures  show  the  dynamic  results  of  these  tests.  Vehicle  attitude  and  power  tracking 
results  are  shown  in  figure  1  (where  1(a)  shows  the  reference  and  attitude  tracking  performance 
of  cases  I,  III,  IV,  and  V,  all  of  which  are  similar,  1(b)  shows  the  attitude  tracking  performance 
of  case  II,  capturing  its  4-channel  quaternion  tracking  error  in  the  4  circles  on  the  plot,  and  1  (c) 
shows  the  reference  power  and  attitude  tracking  case  applicable  to  all  5  cases)  along  with  the 
CMG  singularity  condition  for  all  5  cases  1(d).  It’s  important  to  note  that  all  5  cases  have  the 
same  reference  and  actual  power  profile  plots.  This  means  that  the  cluster’s  unsaturated  wheels 
perfectly  track  the  given  open-loop  power  profile.  This  follows  directly  from  the  cluster’s  VSCMG 
redundancy.  Also,  the  tracking  error  from  case  II  lasts  more  than  200  s.  Unchecked,  this  tracking 
error  can  negatively  impact  precise  imaging  operations. 

These  VSCMG  algorithms  can  be  further  investigated  at  the  actuator  level  by  examining  figure 
2.  This  figure  shows  the  gimbal  angle  and  wheel  speed  responses  of  the  5  cases.  As  is  seen  in  these 
plots,  the  gimbal  and  wheel  responses  for  cases  III,  IV,  and  V  (relating  to  figures  2(c),2(d),2(e),and 
2(f))  are  quite  difficult  to  delineate.  However,  cases  I  and  II,  reflected  in  figure  2(a)  for  the  gimbal 
angles  and  2(b)  for  the  wheel  speeds,  show  differences  in  actuator  response  after  the  key  singu¬ 
larity  (occurs  just  before  40005).  In  case  I,  which  uses  no  singularity  avoidance  algorithm,  the 
gimbal  angles  become  fixed  at  a  singularity  and  do  not  change  for  the  remainder  of  the  maneuver. 
Conversely,  the  wheel  speeds  increase  motion  at  this  singularity  to  account  for  this  locked  config¬ 
uration.  The  result  is  that  the  VSCMGs  become  simple  Momentum  Wheels  for  the  remainder  of 
the  maneuver.  This  result  is  also  undesirable  since  the  wheels  cross  zero  and  begin  to  spin  in  the 
opposite  sense.  In  contrast,  case  II  shows  that  the  GSR  is  effective  at  keeping  the  actuators  away 
from  gimbal  lock  with  the  added  performance  degradation  of  the  aforementioned  attitude  error. 
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4.0  Conclusions 


The  proposed  extended  GSR  and  modified  Schaub  with  null  motion  algorithms  eliminate  small 
attitude  errors  for  all  time  such  that  any  escapable  singularities  occurring  during  a  maneuver  can 
be  mitigated.  In  addition  to  this,  a  given  power  profile  is  simultaneously  tracked  completely 
in  this  case  due  to  cluster  redundancy  based  on  employment  of  the  power  projection  method 
from  [3].  Both  methods  have  a  small  performance  improvement  over  Yoon/Tsiotras’  complex  gra¬ 
dient  method  in  exploiting  each  VSCMG  in  the  cluster’s  CMG-mode,  avoiding  CMG  singularities, 
following  a  desired  3-axis  attitude  trajectory,  and  tracking  a  given  flywheel  battery  power  profile. 
This  result  forms  a  link  between  classical  single-gimbal  CMG  and  modem  VSCMG  theory.  From 
it,  one  can  see  the  relationship  between  the  extended  GSR,  modified  Schaub  VSCMG+null  motion 
steering  law,  and  the  Yoon/Tsiotras  results  from  [5],  each  originally  developed  for  a  different  pur¬ 
pose.  Satellite  designers  can  use  this  result  to  strike  harmony  between  gimbal  and  wheel  motion 
with  important  practical  implications  -  the  ability  to  maintain  precise  pointing  during  agile  slewing 
operations  while  using  VSCMG  flywheels  to  store  (drain)  energy. 
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Appendix  A:  Projected  Power  in  Null  Space  of  Q 


Define  u„  = 


n. 

8 


and  the  operator,  Pn,  which  projects  a  vector,  v  onto  the  nullspace  of 

null 


Q.  Then,  one  can  use  the  step-by-step  equation  set  below  to  derive  the  simultaneous  wheel  and 
gimbal  steering  law,  u.  Notice  that  in  this  equation  set,  represents  any  suitable  generalized 
inverse,  including  Q *  or  Q°. 


Qu  =  Nc 

(23) 

Qun  =  0 

(24) 

Un  =  PnV 
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u  =  Q*NC  +  un 

(27) 

Su  =  S(Q'Nc  +  un)  =  Pc 

(28) 

Sun  —  SPnv  =  Pc-  SQ^Nc 

(29) 

{SPn)v  =  Pc-SQfNc 

(30) 

v  =  (SPn)HPc-SQ'Nc) 

(31) 

v  =  PTnSr{SPnPTnSryX(Pc-S&Nc) 

(32) 

un  =  PnPTnST  (SPnPTnSTYX  {Pc-StfK) 

(33) 

un  =  PnST  (SPnSTyl  (Pc- SQ'NC) 

(34) 

u  =  Q'Nc  +  PnS7  (SPnST)~'  (Pc- SQ^Nc) 

(35) 
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Table  1:  Simulation  Parameters 
Symbol  Value 


Nvc 

4 

0,  deg 

54.75 

co(0),  rad/s 

[0  00] 

P(0) 

[-0.5  0.5  -0.5  0.5] 

8(0),  rad 

r£  .5  .5  2E  1 

U  4  4  4  J 

8(0),  rad/s 

[0  0  0  0] 

n(0),  rpm  [20000  17500  1 5000  1 3000  ] 

I w  j,  kgm2 

diag{0.70,  0.40,  0.40} 

IGj,  kgm2 

diag{0.10  0.10  ,0.10} 

'  15053  3000  -1000  ' 

he,  kgm2 

3000  6510  2000 

-1000  2000  11122 

K,k 

«,  W2jn,Wgjn 


1  0.01  0.02 
0.01  1  0.01 
0.02  0.01  1 
iiag{0.1, 0.1, 0.1}, 0.01 
1,  1000,  10,  0.0001,  le-4 
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Table  2:  Parametric  Comparison  of  Maximum  Attitude  Error,  Momentum  Wheel  Torque,  Wheel 
Speed,  Average  Wheel  Speed,  and  Wheel  Speed  Difference  for  a  Simultaneous,  5000  s  Attitude  Plus 
Power  Tracking  Maneuver 


Case  I 

Case  II 

<jT(P)Pr,  max  att 

3.3015 

3.1430 

4.4373 

4.9126 

804.00 

192.00 

592.00 

182.20 

error,  XI 0-4 
CMG# 

1 

2 

3 

4 

1 

2 

3 

4 

^VwNm 

4.1246 

2.2569 

5.1089 

1.6047 

0.8587 

1.0263 

0.7332 

0.7715 

£W104  RPM 

2.2229 

1.9452 

2.7262 

1.7813 

2.2232 

1.9467 

1.6688 

n,!W,„104  RPM 

1.8201 

1.8209 

%n 

118.86 

184.62 

254.40 

76.699 

25.770 

11.989 

8.3969 

23.982 

Case  III 

<7r(p)Pr,  max  att 

5.6532 

3.9822 

4.421 1 

7.5604 

error,  XI  0~4 
CMG# 

1 

2 

3 

4 

Ww’Nm 

0.8647 

0.7569 

0.6486 

0.5618 

SW,104  RPM 

2.2225 

1.9456 

1.6678 

1.4446 

anom,io4  RPM 

1.8201 

%  Q 

25.768 

11.983 

8.4021 

21.663 

Case  IV 

Case  V 

qT( P)Pr,  max  att 

3.5566 

3.2771 

4.9330 

4.6487 

3.7380 

4.2350 

5.1530 

4.4840 

error,  XI 0”4 
CMG# 

1 

2 

3 

4 

1 

2 

3 

4 

NSmax'  Nm 

0.8640 

0.7560 

0.6480 

0.5616 

0.8574 

0.7525 

0.6427 

0.5575 

amax,\oA  RPM 

2.2230 

1.9451 

1.6672 

1.4449 

2.2234 

1.9452 

1.6671 

1.4446 

<W104RPM 

1.8201 

1.8201 

%a 

25.768 

11.983 

8.397 

20.622 

26.212 

7.820 

11.295 

24.304 
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a)  Attitude  Tracking  for  Cases  I,  III,  IV,  V 
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b)  Attitude  Tracking  for  Case  II 


c)  Reference  and  Actual  Power  Profile  for  Cases  I,  III,  IV,  V  d)  Singularity  Condition  for  Case  I-V 

Figure  1:  Attitude,  Power  Tracking,  and  Singularity  Condition 
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e)  Case  III  (Extended  GSR)  Gimbal  Angles 
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f)  Case  III  (Extended  GSR)  Wheel  Speeds 


Figure  2:  Actuator  States  Comparison 
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